Physics of Psychophysics: Stevens and Weber-Fechner laws are transfer functions of 
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Sensory arrays made of coupled excitable elements can improve both their input sensitivity and 
dynamic range due to collective non-linear wave properties. This mechanism is studied in a neural 
network of electrically coupled (e.g. via gap junctions) elements subject to a Poisson signal process. 
The network response interpolates between a Weber-Fechner logarithmic law and a Stevens power 
law depending on the relative refractory period of the cell. Therefore, these non-linear transforma- 
tions of the input level could be performed in the sensory periphery simply due to a basic property: 
the transfer function of excitable media. 

PACS numbers: 05.45.Ra, 05.45.Xt, 87.10.+e, 87.18. Sn 
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A very common trade-off problem found in the biology 
of sensory mechanisms (and sensor devices in general) is 
the competition between two desirable goals: high sensi- 
tivity (the system ideally should be able to detect even 
single signal events) and large dynamic range (the system 
should not saturate over various orders of magnitude of 
input intensity). In physiology, for example, broad dy- 
namic ranges are related to well known psychophysical 
laws pi 0|: the response R of the sensory system may 
be proportional not to the input level I but to its loga- 
rithm, R oc In/ (Weber-Fechner Law) or to a power of 
it, R oc 7 a , (a < 1) (Stevens Law). 

Most of the attempts to explain these psychophysics 
laws consist basically in top-down approaches trying to 
show that they could be derived from some optimiza- 
tion criterium for information processing [3- ■ In this 
work we use a botton-up, statistical mechanics approach, 
showing how these laws emerge from a microscopic level. 
Indeed, they are generic transfer functions of excitable 
media subjected to external (Poisson) input. Of course, 
this does not explain "why" these laws have been adopted 
by Biology (some optimization criterium may be relevant 
here), but explains why Biology uses excitable media to 
implement them. 

Receptor cells of sensory systems are electrically cou- 
pled via gap junctions 0, 13 • However, the functional 
roles of this electrical coupling are largely unknown. Here 
we report a simple mechanism that could increase at 
the same time the sensitivity and the dynamic range of 
a sensory epithelium by using only this electrical cou- 



pling. The resulting effect is to transform the individual 
linear-saturating curves of individual cells into a collec- 
tive Weber-Fechner like logarithmic response curve with 
high sensitivity to single events and large dynamic range. 
We also observe a change to power law behavior (Stevens 
Law) if relative refractory periods are introduced in the 
model. 

Although the phenomenon discusssed in this work 
could be illustrated at different modeling levels we 
have chosen here to work with the simplest elements: 
cellular automata (CA). The simplicity of the model sup- 
ports our case that the mechanisms underlying the de- 
scribed phenomena are very general. To confirm this 
picture, we also present preliminary results for neurons 
modeled by the Hodgkin-Huxley equations. 

The n-state CA model is an excitable element contain- 
ing two ingredients: 1) a cell spikes only if stimulated 
while in its resting state and 2) after a spike, a refractory 
period takes place, during which no further spikes occur, 
until the cell returns to its resting state. Denoting the 
state of the i-th cell at time t by Xi(t) £ {0, 1, . . . , n— 1}, 
the dynamics of the proposed CA can be simply described 
by the following rules: 



0, then Xi(t + 1) = hi(t), where hi € 
2. If Xi(t) ^ 0, then x t (t + 1) = [xi(t) + 1] mod n. 



1. If Xi (t) 
{0,1}. 
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Interpretation of the above rules is straightforward: a cell 
only responds to stimuli in its resting state (xi = 0). If 
there is no stimulus (hi — 0), it remains unchanged. In 
case of stimulus (hi = 1), it responds by spiking (xi = 1) 
and then remaining insensitive to further stimuli during 
n — 2 time steps (xi £ {2, . . . , n — 1}). 

In what follows, we assume that the external input 
signal Ii(t) arriving on cell i at time t is modeled by a 
Poisson process of supra-threshold events of stereotyped 
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FIG. 1: Time evolution for n — 5: (a) V(x(t)) for a single 
uncoupled cell (solid lines) and stimuli h(t) (bars) at r = 
100 events/sec with Vo = 5 and k = 0.3; (b) system with 
L — 50 coupled cells at r = 10 events/sec: Xi — 1 (filled 
circles), 2 < Xi < n — 1 (open circles) and Xi = (white 
background). Arrows indicate events to be considered in more 
detail subsequently. 



unit amplitude: 7j(t) = S [t, t„ J where S(a,b) is 

(i) (i) 

the Kronecker delta and the time intervals t n+1 — t n 
are distributed exponentially with average (input rate) 
r, measured in events per second. For uncoupled cells, 
we have then simply hi(t) = 6 (Ii(t), 1)- 

In order to visualize the effect of the refractory period, 
we mimick the behavior of the spike of a neuron by map- 
ping the automaton state into an action potential wave 
form 

V(xi) = V {8(xi,l)-[l-5(xi,0)][l-5{xi,l)] 

Notice that V plays no role whatsoever in the dynamics. 
Fig. da) shows the behavior of V(xi(t)) for an uncou- 
pled 5-state automaton. We observe that stimuli that 
fall within the refractory period go undetected, and in 
the absence of stimuli the automaton eventually returns 
to and stays at its quiescent state x\ = 0. Since a typical 
spike lasts the order of 1 ms, this provides a natural time 
scale of 1 ms per time step, which will be used throughout 
this paper. 



Response of uncoupled receptor cells is shown in Fig. 
(thick lines on top panels). We draw input signals at rate 
r per cell and measure the average firing rate / (spikes 
per second per cell) of the n-state automata over a suf- 
ficiently long time. In the low rate regime the activity 
of the uncoupled cells is proportional to the signal rate. 
If the rate increases, there is a deviation from the linear 
behavior due to the cell's refractory time A„ = n x 10~ 3 
seconds. The single-cell response / is extremely well fit- 
ted by a linear-saturating curve /„ [Fig |2f a) and (b)] : 

/ n (r)=r/(l + rA„), (2) 

which can be deduced from the fact that the firing rate is 
proportional to the rate discounting the refractory inter- 
vals, /„ = r(l — /„A„). The same result can be obtained 
by a stationary mean field solution of the uncoupled cells. 

How to improve the sensitivity for very low rates? If 
we consider the response R (spikes per second) of the to- 
tal pool of L independent cells, we have R = Lf w Lr, 
so increasing L increases the total sensitivity of the ep- 
ithelium. Although certainly useful, this scaling is trivial 
since the efficiency of each cell remains the same. 

Coupled excitable cells (say, via gap junctions) are an 
example of excitable media that supports the propaga- 
tion of nonlinear waves • Here we show that the forma- 
tion and annihilation of these waves enhance the sensi- 
tivity and, at the same time, extends the dynamic range 
of a sensory epithelium. We couple L cellular automata 
in a chain by defining the local input as 

h i (t) = l-[l-S(I i (t),l)} [] [l-6(x i+j (t),l)} , (3) 

3=±1 

i.e. hi(t) will be nonzero whenever either of i's neighbors 
are spiking and/or the external input is nonzero. This 
kind of coupling models electric gap junctions instead of 
chemical synapses because it is fast and bidirectional. 

A sample of the resulting chain dynamics is shown in 
Fig. ^b). Due to coupling, single input events create 
waves that propagate along the chain, leaving behind a 
trail of refractoriness (of width n — 1) which prevents new 
spikes from reappearing immediately. More importantly, 
refractoriness is responsible for wave annihilation: when 
two wave fronts meet at site i they get trapped because 
the neighboring sites have just been visited and are still in 
their refractory period. This is a well known phenomenon 
in excitable media and occurs in the CA model Mn > 3. 
Notice that the overall shape of two consecutive wave- 
fronts are correlated (see Fig. 1), denoting some kind 
of memory effect, a phenomenon observed previously by 
Chialvo et al. Q and Lewis and Rinzel . 

Due to a chain-reaction mechanism, the spike of a sin- 
gle receptor cell is able to excite all the other cells. The 
sensitivity per neuron has thus increased by a factor of 
L. This can be clearly seen in Fig. [21 which shows the 
average firing rate per cell F in the coupled system (top 
panels), as well as the amplification factor A = F j f (bot- 
tom panels). This is a somewhat expected effect of the 
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FIG. 2: Firing rates / and F (top) and amplification factor 
A (bottom) vs. input rate r for n = 3 and varying L [(a) and 
(c)] and for L — 5000 and varying n [(b) and (d)]. Thick lines 
(top panels) show f n {r) as in Eq.|3 
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FIG. 3: F x A vs. input rate r for L — 5000 (open symbols) 
and L — 200 (filled symbols) for different values of n. A 
L = 25000 curve for n — 50 (crosses) shows no difference to 
the L — 5000 case. Straight lines are intended as a guide to 
the eye. Inset: F(r) for the Hodgkin-Huxley system. 



coupling: neuron j is excited by signal events that arrive 
not only at neuron j but elsewhere in the network. 

More surprising is the fact that the dynamic range (the 
interval of rates where the neuron produces appreciable 
but still non-saturating response) also increases dramat- 
ically. This occurs due to a second effect, which we call 
the self-limited amplification effect. Remember that a 
single spike of some neuron produces a total of L neu- 
ronal responses. This is valid for small rates, where in- 
puts are isolated in time from each other. However, for 
higher signal rates, a new event occurs at neuron k be- 
fore the wave produced by neuron j has disappeared. If 
the initiation site k is inside the fronts of the previous 
wave [e.g. the events signaled by arrows in Fig. ^b)], 
then two events produce 2L responses as before. But if 
k is situated outside the fronts of the j-initiated wave [as 
in the first input events shown in Fig. Hfb)], one of its 
fronts will run toward the j-wave and both fronts will 
annihilate. 

Thus, two events in the array have produced only L 
excitations (that is, an average of L/2 per input event). 
So, in this case, the efficiency for two consecutive events 
(within a window defined by the wave velocity and the 
size L of the array) has been decreased by half. If more 
events (say, m) arrive during a time window, many fronts 
coexist but the average amplification of these m events 
(how many neurons each event excites) is only of order 
L/m. 

Therefore, although the amplification for small rates 
is very high, saturation is avoided due to the fact that 
the amplification factor decreases with the rate in a self- 
organized non-linear way. The amplification factor A 
shown in Figs. [2Jc) and |2Id) decreases in a sigmoidal 
way from A = O(L) for very small rates (since a single 
event produces a global wave) to A = 1 for large rates, 
where each cell responds as if isolated since waves have 
no time to be created or propagate. 

The role of the system size L for low input rates be- 



comes clear in Fig.EJc): the larger the system, the lower 
the rate r has to be in order for the amplification factor 
to saturate at O(L). In other words, we can think of a 
decreasing crossover value r% (L) such that the response is 
well approximated by F(r) = Lf(r) ks Lr for r <C ri(L). 
In this linear regime consecutive events essentially do not 
interact. Larger system sizes increase not only the overall 
rate of wave creation (~ 1 — (1 — r) L ) but also the time it 
takes for a wave to reach the borders and disappear. In 
the opposite limit of large input rates, the behavior of the 
response is controlled by the absolute refractory period 
A, as shown in Fig. EI F and / saturate at r2 = I/A for 
/ > Tii independently of the system size. 

So what happens for intermediate input rates, i.e. 
r i ^ r ^ r 2? The answer is a slow, Weber- Fechner- like 
increase in the response F, as can be seen in Fig. 01 The 
logarithmic dependence on r is a good fit of the curves 
for about three decades. 

Motivated by results obtained with more realistic el- 
ements Q we introduced a relative refractory period in 
our CA model. We first define a time window M after a 
spike during which no further spikes can occur (absolute 
refractory period) . In the following n — M — 2 steps (rela- 
tive refractory period), a single input does not produce a 
spike but two or more inputs can elicit a cell spike if they 
arrive within a temporal summation window r (details of 
this model will be described in a forthcoming full paper). 
This ingredient produced the appearance of a power law 
F(r) curve (Stevens Law PJ0)) as shown in Fig. 4. No- 
tice that the exponent depends on the relative refractory 
period. The appearence of a power law transfer function 
is a robust effect also observed in coupled maps systems 

We may confirm the generic character of the self- 
regulated amplification phenomenon by performing sim- 
ulations using biophysically detailed cell models, for ex- 
ample a network of Hodgkin-Huxley (HH) elements with 
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FIG. 4: Neuronal "Stevens Law" F oc r a in automata which 
takes temporal summation effects into account (see text for 
details) . Firing rate F vs. input rate r for a CA with n states 
and an absolute refractory period of M = 3 time steps. Filled 
circles: n = 15, r = 10, a = 0.38; open circles: n = 100, 
r = 80, a = 0.44. 
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FIG. 5: (a) Firing rate for coupled (F, filled circles) and un- 
coupled (/, open circles) systems and (b) amplification factor 
A = F/f vs. input rate r for Hodgkin-Hukley neurons for 
L = 200. 



the standard set of parameters given in 0] connected 
via gap junctions of 30A/f2. Preliminary results show 
that this system exbits the same qualitative behavior of 



the simple CA model (see inset of Fig. |3] and Fig. 0. 
More detailed results will be reported elsewhere. 

Concerning the functional role of gap junctions for sig- 
nal processing, it has been recognized that they pro- 
vide faster communication between cells than chemical 
synapses and play a role in the synchronization of cell 
populations [ljj . Here we are proposing another func- 
tional role for gap junctions: the enhancement of the 
dynamic range of neural networks. 

There is considerable debate about what is the most 
appropriate functional law to describe psycophysical re- 
sponse: Weber-Fechner, Stevens or some interpolation 
between the two Q- Our results suggest that proper- 
ties of excitable media could be a bottom-up mechanism 
which can generate both laws, and a cross-over between 
them, depending on the presence of secondary factors like 
the relative refractory periods and temporal summation. 

We can even make two more specific predictions which 
are easily testable experimentally: 1. The larger the rel- 
ative refractory period (e.g., due to slower hyperpolar- 
izing currents) of sensory epithelia neurons, the larger 
the exponent of Stevens Law; 2. For sufficiently low in- 
put rates, the sensory epithelium response will be always 
linear (a = 1). 

This mechanism for amplified but self-limited response 
due to wave annihilation promotes signal compression, is 
a basic property of excitable media and is not restricted 
to one dimensional systems. We conjecture that the same 
mechanism could be implemented at different biologi- 
cal levels, from hippocampal networks (where axo-axonal 
gap junctions have been recently reported 11] and mod- 
eled by a CA similar to ours) to excitable dendritic 
trees in single neurons 0, 0] . This signal compression 
mechanism could also be implemented in artificial sensors 
based on excitable media. 
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